## Requires one argument: .igv file

import sys

TA = 0; Reads = 0; HitSites = 0

for line in open(sys.argv[1]):
    split=line.split()
    if len(split) > 3:
        if split[3][0].isdigit() == True:
            TA += 1
            Reads += float(split[3])
            if float(split[3]) > 0: HitSites += 1
    
print 'TA = %d' % (TA)
print 'Reads = %d' % (Reads)
print 'Sites Hit = %d' % (HitSites)
